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The time evolution of a single particle in a harmonic trap with time dependent frequency u(t) 
is well studied. Nevertheless here we show that, when the harmonic trap is opened (or closed) as 
function of time while keeping the adiabatic parameter fi = [duj(t)/dt]/uj 2 (t) fixed, a sharp transition 
from an oscillatory to a monotonic exponential dynamics occurs at (X = 2. At this transition point 
the time evolution has a third-order exceptional point (EP) at all instants. This situation, where an 
EP of a time-dependent Hermitian Hamiltonian is obtained at any given time, is very different from 
other known cases. Our finding is relevant to the dynamics of a single ion in a magnetic, optical, or 
r/trap, and of diluted gases of ultracold atoms in optical traps. 



Exceptional points (EP) are degeneracies of non- 
Hermitian Hamiltonians [IJ [2], associated with the co- 
alescence of two or more eigenstates. The studies of 
EPs have substantially grown since the pioneering works 
of Carl Bender and his co-workers on 'PT-symmctric 
Hamiltonians [3]. These Hamiltonians have a real spec- 
trum, which becomes complex at the EP. However, VT- 
symmetry is not required to obtain an EP point, as in 
the case of a coalescence between two resonant states, 
leading to self-orthogonal states jlHS]. 

The physical effects of second-order EPs have al- 
ready been demonstrated in different types of experi- 
ments. See for example the effect of EPs on cold atoms 
experiments [7 1, on the cross sections of electron scatter- 
ing from hydrogen molecules [8], and on the linewidth of 
unstable lasers [9]. More direct realizations of EPs in mi- 
crowave experiments are given in Ref. [101 lllj and in op- 
tical experiments in Ref. [T2] . For theoretical studies that 
are relevant to these experiments see for example [HI 1131 - 
fTB"] . In addition, theoretical studies predict significant 
effects of second-order EPs on the photoionization of 
atoms [T5H2I] and the photodissociation of molecules [221- 
[2J]. 

The possibility of higher-order EPs (where more than 
two eigenstates of the non-Hermitian Hamiltonian co- 
alesce at the EP) has been discussed in the literature 
for time independent VT symmetric Hamiltonians (see 
Ref. |25l [25] and references therein). The main effect of 
EPs (of any order) on the dynamics of "PT-symmetric 
systems is the sudden transition from a real spectrum to 
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a complex energy spectrum associated with gain and loss 
processes. [H]. 

All above mentioned studies on the effects of EPs are 
related to non-Hermitian time-independent Hamiltoni- 
ans. Note that non-Hermitian Hamiltonians can be ob- 
tained from Hermitian Hamiltonians by imposing outgo- 
ing boundary conditions on the eigenfucntions or includ- 
ing complex absorbing potentials [5]. This approach al- 
lows the description of resonance phenomena in systems 
with finite-lifetime metastable states. 

Other studies considered time-periodic Hamiltonians 
where the EPs are associated with the quasi-energies 
of the Floquet operator which can be represented by a 
time-independent non-Hermitian matrix (see for exam- 
ple one of the first studies of EP in atomic physics in 
Ref.|19j). EPs were also studied in non-periodic sys- 
tems in the context of Landau-Zener-Majorana transi- 
tions, where the EP was obtained only after analytic con- 
tinuation of the actual Hamiltonian [27j. Finally, time- 
dependent EPs have been used to control the quantum 
evolution of non-Hermitian systems [28 . 

In this paper we show that an EP for the time evo- 
lution of a Hermitian time dependent Hamiltonian can 
be obtained. Our method relies on the re-scaling of the 
time-axis, allowing to map the time-dependent problem 
to an effective time-independent one, with non-unitary 
evolution. 

The harmonic oscillator system with changing fre- 
quency in the Heisenberg picture - The ID harmonic os- 
cillator system with changing frequency is defined as 

H=^-p 2 + \muj 2 {t)x\ (1) 
2m 2 

where m is the mass of the particle, and p and x 
are respectively the momentum and position operators. 
Despite its simplicity, this time-dependent Hermitian 
Hamiltonian can associated with a third-order excep- 
tional point. 

We study the model ([I]) in the framework of Ref. [291 
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[30], where it is shown that, due to the closed commuta- 
tion relations between the operators p 2 ,x 2 ,px + xp, the 
model forms a SU(1,1) algebra. As a basis set for this 
algebra we choose the Hamiltonian, the Lagrangian and 
the x — p anti-commutator: 



01 = H 

2 = L = H — muj 2 (t)x 2 

3 = D = uj(t)(xp + px)/2 



(2) 
(3) 
(4) 



Any commutator between operators in the algebra can 
be expressed as a linear combination of these operators: 
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(h.O, , 1,2,3 }C',,P 



k=l 



(5) 



where the Cj fe are the structure factors of the SU(1,1) 
algebra [5Tj. The Heisenberg picture for the dynamics 
which is associated with the operators Oj is described as 



where j = 1, 2, 3. 

These equations are explicitly given by 



(6) 



at 



dH 
dt 

d f i_ 



nuiujx 2 = cufi(H — L) 
[H, £] + # = -2uD - LJfj,(H - L) 



dt,' 

jr[H,D] + ^ = 2uL + unD 
Here we defined the dimensionless "adiabatic parameter" 



(7) 



/' 



W 3(t) 



doj 
~dl' 



([7| conserve the "Casimir" 
-L 2 (t)-D 2 (t)}/uo 2 {t) by sat- 



The equations of motion 
operatorp] 6{t) = [H 2 {t)- 
isfying dC '/dt — 0. 

In what follows we will focus on the specific case of 
\i = const, corresponding to the frequency profile 



w(t) 



"(0) 
l-/iw(0)i' 



(9) 



In experiments the harmonic trap is varied between two 



extreme values, uj open and oj c i 



The compression fac- 



For positive values of the 

open and u(tf) = UJrfosed- 



open 



. In 



tor is given by uj c i ose d 
adiabatic parameter w(0) = cu, 
For negative values w(0) = uj c i ose( i and ui(tf) 
both cases t f = \n\(uj^ p \ n - u^J). 

The parameter [i sets the degree of adiabaticity of the 
process. For fi — > 0, the dynamics is perfectly adiabatic 
and the system follows the eigenvalues of the instanta- 
neous Hamiltonian. In contrast, for [i — > ±00, the change 
of the Hamiltonian is so fast that the system does not 
have time to change at all. As we will show, these two 
limits are separated by an exceptional point. A similar 



effect is known to occur in the vicinity of quantum crit- 
ical points (See for example Ref . [33] ) and is here shown 
in time-dependent non-periodic harmonic traps. 

It is convenient to introduce the dimensionless time 
variable r = log(u;(£)/cj(0)), satisfying dr = 

u>(t)dt, and rewrite ^ as 



dO(-r) 



= ^HoiaO(r)= (zMl + ^Sd 



0(r) (10) 



where O = {H, L, _D}-components, I is the unit matrix 
and the traceless operator H]^[® is defined by 



n 



Tr=0 
Hcis 




(11) 



We can further simplify Eq.( 10 1 by performing the trans- 
formation 



0(r(t)) 



dO(r) dO(r) 



u(t) 



(r(t)) 

dr dr 

(12) 

The resulting equation of motion idO/dr = H]^~®0 
is equivalent to a time dependent Schrodinger equa- 
tion with a non-Hcrmitian time-independent Hamilto- 
nian. The matrix T-L]^[ s is 7 ? 7~-svmmetric pM] and its 
three eigenvalues 



(13) 



(8) ( 



are real for < 2 

vectors are given by, 



The corresponding eigen- 
v = (1,0,-^/2), v± = 

A?, -2) //*. 

In contrast to the Schrodinger equation, the population 
of the eigenvectors in a physical state is not arbitrary, 
but must satisfy several constraints. For example, for 
I /1 1 < 2, the eigenvectors v + and u_ are complex and 
any physical state must populate them with an equal 
weight, in order to keep the expectation values H, D, 
and L real. In addition, v + and u_ have a zero Casimir 
constant (C) = H 2 — L 2 —D 2 = 0. Due to the uncertainty 
relation (C) > h 2 /A any physical state must necessarily 
populate the eigenstate vq with non-zero weight as well. 
Thus, for a generic initial state, we expect more than 
one eigenvector to be occupied, leading to an oscillatory 
behavior that we describe below. 

The third-order EP for the time evolution operator of 
the Hermitian time- dependent harmonic oscillator - The 
matrix H^^ has a third order EP at = 2. At this 
point all three eigenvalues and the corresponding eigen- 
vectors (in the HLD space) coalesce. As a consequence, 



u 



Tr=0 



(M = ±2) 



while 



n 



Tr=0 



(m = ±2; 



0, 



demonstrating that the present EP is of third order. 

At the EP the matrix "H^^ has one single eigenvec- 
tor, which is "self-orthogonal" . To show this property, 
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it is necessary to multiply the right and left eigenvectors 
of the non-symmetric "Hh^ - At the EP the right eigen- 
vector is (1,0,-1) while the left eigenvector is (1,0,1). 
Their product is equal to zero showing that the eigen- 
vector is "self-orthogonal" (see for example Chapter 9 in 
Ref.0). 

Let us first discuss the situation where the initial state 
is Gaussian (e.g., a ground state of the harmonic oscilla- 
tor trap at t = or a thermal state). The evolution in 
time is described by a Gaussian Wigner distribution with 
elliptic contour plots. Using the normalized coordinates, 
x(moj) 1 / 2 and p/imio) 1 / 2 , the standard deviation of the 
narrow axis (N) and the wide axis (W), are given by: 



H±y/L 2 + D 2 



W.N 



(14) 



Note that by virtue of the Casimir constant a most gen- 
eral state satisfies <jn<Jw > 1/2- For uj = const (p = 0) 
the evolution simply mixes L and D leaving \J L 2 + D 2 
constant. Hence even though for u> — const the distri- 
bution rotates in phase space and changes the variance 
of position and momentum, the width of the narrow and 
wide axis of the distribution remains fixed. 

For < 2 time oscillations of measurable quantities 
are obtained. A convenient measure to capture this os- 
cillatory dynamics is given by the ratio: 



<jw_ 



(15) 



where by definition p > 1. The time evolution of p is 
shown in Fig. la as function of the time-dependent com- 
pression factor uj(t)/uj(0) in a log scale. Since the loga- 
rithm of the compression factor is equal to the new time 
variable t (multiplied by p), in this scale p(t) shows pe- 
riodic oscillations. In the original time variable t, the 
system displays non-periodic oscillations with the same 
amplitude. The periodicity in the r time coordinate is 
determined by the eigenvalue difference: 



T T = 



2tt 



2tt 



\(E+-EJ) JI- 



(16) 



The visibility of the fringes pattern that appear in Fig. 
la is given by the simple expression: 



V = 



Pn 



Pmin | | 

Pmin 2 



(17) 



At \p\ = 2 the visibility becomes one, and the oscillations 
disappear for \p\ > 2. 

We now comment on the time evolution starting from 
a generic non-Gaussian initial state. Of course, in this 
case H,L and D are not sufficient to completely deter- 
mine the state. Nevertheless for p — const the expec- 
tation values of H,L and D follow exactly the dynam- 
ics described above. It is possible to show that the 



pD(0) + L(0) 2 /(4ff(0)) > 0. In particular, this condi- 
tion is fulfilled by any initial state that is stationary with 
respect to the initial Hamiltonian H(0). 

A similar transition between oscillatory and 
monotonous behavior appears in other properly scaled 
quantities (see Eq 12), like (H — L)/uj{t) = muj(t) (x 2 ) 
and {H + L) /w{t) = (p 2 ) /(muj(t)), as shown in Fig. lb. 
However note that for initial states which are not 
eigenstates of the harmonic trap these quantities display 
additional trivial oscillations that show up even for 
lj = const. 

The EP corresponds to a transition from an under- 
damped to an over-damped harmonic oscillator- We now 
present a different approach which will clarify the relation 
between the present problem and energy-dissipative sys- 
tems. Our approach is based on the equivalence between 
quantum and classical evolution of quadratic Hamiltoni- 
ans. To reproduce the quantum mechanical results one 
simply needs to complement the classical equations of 
motion by stochastic initial distributions, given by the 
Wigner transform of the initial state. 

In our case, the relevant equation of motion is New- 
ton's law 



dt 2 



■u 2 {t) 



x(t) = 



(18) 



By applying the transformation dr — oj(t)dt, or 4? = 
w(i) J^, we obtain 

a¥ x = jt {Lj{t) i x) = L0 ' {t) i f+ ^ {t) i x (19) 

In the specific case p = const, the equation of motion 
then becomes 



dr 2 ^ dr 



x(t) = 



(20) 



Here we obtain the time-independent non-Hermitian 
equation of motion of a damped harmonic oscillator. 
The EP \p\ =2 corresponds to the transition between 
an under-damped and over-damped oscillator, as can be 



seen by the Fourier transform of ( 20 1 , leading to 



A, 



ip±y/4- p 2 



(21) 



expression (17) is valid for any initial state satisfying 



For \p\ > 2 both eigenfrcquencies are pure imaginary, 
leading to the disappearance of the oscillatory behavior. 

Physical realization of the third-order EP in experi- 
ments - Although our Hamiltonian can be realized in any 
controllable Harmonics trap (optics, plasma, ...), we will 
consider here the case of either a single particle (ion|35j). 
or a dilute atomic cloud 36 -38 , in time-dependent con- 
fining traps. The realization with atomic clouds allows 
the measurement of expectation values in a single-shot 
experiment. Complications of the dynamics due to the 
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FIG. 1: (a) The ratio, p, between the narrow and wide axis 
in the Gaussian Wigner distribution as a function of the time 
dependent compression factor w(t)/w(0) for different values 
of the adiabatic parameter /i > 0. For /j, < the same 
plot is obtained where now the compression factor is taken 
as Lj(0)/u)(t). At the exceptional point (EP) u = 2 the dy- 
namics changes from oscillatory to monotonous, (b) Same 
plot for the variance of the position operator normalized by 
the instantaneous frequency, 2{x 2 (t))uj(t). The initial state 
is the ground state of the Hamiltonian Q at t — where 
(a; 2 ) = l/(2mw(0)). 



atom-atom interactions can be avoided (minimized) by 
setting the atomic scattering length to zero in the vicin- 
ity of a Feshbach resonance |39j. 

The procedure to observe the EP effect on the dy- 
namics of time-dependent Hermitian Hamiltonian is as 
follows: 

(1) Equilibrate the matter in a harmonic trap character- 
ized by the frequency w(0). 

(2) Vary the frequency of the trap as function of time 
from w(0) to w(i/), while imposing a constant adiabatic 
parameter /z, as shown in Eq|9] 

(3) At time {0 < t n < i/}n=i,2...,Af measure either the 
spatial distribution of matter inside the trap[38 l I40H44"] . 
or the momentum distribution in a time of flight experi- 
ment, by suddenly turning off the trap. 

(4) Scale the variance of the measured data (posi- 
tion or momentum) by the instantaneous frequency: 
(x 2 (t)) -> u(t)(x 2 (t)) and (p 2 (t)) -> {x 2 (t))/u{t). 
These quantities are then plotted as function of the 
compression factor u(t)/uj(0) for different values of fx as 
in Fig. lb, showing oscillatory behavior for any |/x| < 2, 
and exponential behavior for > 2. 



To avoid the need to use an eigenstate of the Hamil- 
tonian at t = 0, one should plot the ratio between the 
narrow and wide axis of the Wigner distribution of the 
propagated wavepacket on and ow, defined above. This 
ratio can be measured as follows. At time i„, rather 
than performing a direct measurement, we propose to 
keep the frequency of the trap unchanged at uj = oj{t n ) 
and to measure the variance of the position (or momen- 
tum) as function of time. The minimum and the maxi- 
mum of mu](t n )(x 2 (t)} are respectively a 2 N and o w . This 
method is perhaps more time-consuming, but guarantees 
the independence of the result on the initial preparation 
provided that ^D(O) + L(0) 2 /(AH (0)) > 0. Note that an 
initial stationary state (either an eigenstate of the initial 
Hamiltonian or a thermal state) satisfies this condition. 

Concluding remarks - The dramatic effect of EPs of 
non-Hcrmitian time-independent Hamiltonian systems 
on the dynamics is in the focus of recent theoretical 
and experimental studies in various fields of physics 
(as for example in optical or microwave experiments 
where the material has a complex index of refraction). 
Here we show that the dynamics of a system described 
by a time- dependent Hermitian Hamiltonian can be 
strongly affected by a third-order EP of an effective time- 
independent Hamiltonian. The fact that the dynamics of 
the Hermitian time-dependent harmonic oscillator can be 
explained by the existence of an EP at all instants shows 
the richness of the dynamics of one of the most basic 
model Hamiltonians, which constitutes a milestone in a 
large variety of fields in physics. Our finding is both in- 
teresting for fundamental theoretical reasons and useful 
to control the dynamics of a single-ion and of a diluted 
BEC in a time-dependent traps. 
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daily basis among the permanent and visiting members 
of ITAMP. 
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